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Abstract 



The observation that the existence of the amazing reality and discreteness 
of the spectrum need not necessarily be attributed to the Hermiticity of the 
Hamiltonian is reemphasized in the context of the non-Hermitian Dirac and 
Klein-Gordon Hamiltonians. Complex Coulombic potentials are considered. 
In one of the first explicit studies of the non-Hermitian Schrodinger Hamiltonians, Caliceti 
et al [1] have considered the imaginary cubic oscillator problem in the context of perturbation 
theory. They have offered the first rigorous explanation why the spectrum in such a model 
may be real and discrete. Only many years later, after being quoted as just a mathematical 
curiosity [2] in the literature, the possible physical relevance of this result reemerged and 
emphasized [3]. Initiating thereafter an extensive discussion which resulted in the proposal 
of the so called VT — symmetric quantum mechanics by Bender and Boettcher [4]. 

The spiritual wisdom of the new formalism lies in the observation that the existence of 
the real spectrum need not necessarily be attributed to the Hermiticity of the Hamiltonian. 
This observation has offered a sufficiently strong motivation for the continued interest in 
the complex, non-Hermitian, cubic model which may be understood as a characteristic 
representation of a very broad class of the so-called pseudo-Hermitian models with real 
spectra. 

In such non-Hermitian settings, new intensive studies employed , for example, the idea 
of the strong coupling expansion [5], the complex WKB [6], Hill determinants and Fourier 
transformation [7], functional analysis [8], variational and truncation techniques [9], linear 
programing [10], pseudo-perturbation technique [11,12], ..etc (cf [13-15]). However such 
studies remain in the context of Schrodinger Hamiltonian and need to be complemented by 
the non-Hermitian setting of Dirac and Klein-Gordon Hamiltonians. Starting, say, with our 
forthcoming oversimplified generalized complex Coulombic examples. 

A priori, a generalized Dirac - Coulomb equation for a mixed potential consists of a 
Lorentz-scalar Coulomb-like and a Lorentz-vector Coulomb potentials. Whilst the former is 
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added to the mass term of the Dirac equation, the minimal coupling is used, as usual, for 
the latter. The ordinary ( Hermitian) Dirac Hamiltonian is exactly solvable in this case ( 
cf. e.g., [16,17]). In fact the exact solution to Dirac equation for an electron in a Coulomb 
field was first obtained by Darwin [18] and Gordon [19]. 

The key idea is that instead of solving Dirac-Coulomb equation directly, one can solve the 
second-order Dirac equation [16-22] which is obtained by multiplying the original equation, 
from the left, by a differential operator. The second-order equation is similar to Klein- 
Gordon equation in a Coulomb field. The latter reduces to a form nearly identical to that of 
the Schrodinger equation and its solution can thus be inferred from the known non-relativistic 
solution. 

In what follows we recycle the modified similarity transformation ( used by Mustafa et al 
[17] ) and obtain exact solutions for the non- Hermitian generalized Dirac and Klein-Gordon 
Coulomb Hamiltonians. Although this problem might be seen as oversimplified, it offers a 
benchmark for the yet to be adequately explored non- Hermitian relativistic Hamiltonians. 

For a mixed scalar and electrostatic complex Coulombic potentials, i.e. m — > m — iA 2 /r 
and V(r) = —iAx/r, the Dirac Hamiltonian reads ( with the units k — c — 1) 



where the Dirac matrices a and (3 have their usual meanings. With the similarity transfor- 
mation 



H = a ■ p + (5{m — iA 2 /r) — iAi/r, 



(1) 



S = a + ib(3a-r ; S l = 



a — ib (3 a ■ f 



(2) 



a 2 -b 2 



applied to Dirac equation one gets 



H 



' = SHS 



-l 



(3) 
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where f is the unit vector f/r and a and b are constants to be determined below. For 
the above central problem, the transformed wave function is given by 



i R(r) ¥ jm 
Q(r)a-r& jm 



(4) 



In a straightforward manner one obtains, through E^>' = S H S two coupled equations 
for R(r) ( the upper component) and Q(r) ( the lower component): 



IK % Ai 
\d r H 1 cosh 9 H sinhO + Esinh9] R(r) = £i(r) Q(r), 

ry ry ry 



(5) 



with 



IK iAi 

\d r + cosh 6 sinh 6- E sinh 6]Q(r) = £ 2 (r) R(r), 

r r r 



iAi i A\ K 

£i(r) = m 1 cosh 6 H sinh 9 + E cosh 9, 

r r r 



(6) 



(7) 



£ 2 (V) — m — ■ — - — - — -cosh 9 sinh 9 — E cosh 9. 

r r r 



(8) 



Where K = u (j + 1/2), Q = ^1 for I = j + Co/2, cosh9 = (a 2 + b 2 )/(a 2 - b 2 ), and 
sinh 9 = 2ab/{a 2 -b 2 ). 

Incorporating the regular asymptotic behaviour of the radial functions near the origin; 
i.e. R(r) — > air 7-1 and Q(r) — > a 2 ^ 7 1 as r — > 0, and neglecting all constant terms propor- 
tional to mass and energy, one obtains 



7 = JK* + A*-A*. 



(9) 
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The negative sign of the square root has to be discarded to avoid divergence of the wave 
functions at the origin. 

It is obvious that one has the freedom to proceed either with the upper radial compo- 
nent R(r) or with the lower component Q(r). We shall, hereinafter, work with the upper 
component and determine sink 9 and cosh 9 ( hence the constants a and b) by requiring 



-i A 2 + i A\ cosh 9 + K sinh 9 = 0, 



(10) 



K cosh 9 + i Ai sinh 9 = cu 7, 



(11) 



This requirement yields 



sinh 9 = — iuj 



„ [Ax 7- \K\A 2 



[K 2 + A 2 ] 

Equations (5) and (6) would, as a result, imply 



, cosh 9 



1^17 + ^1^2 

[K 2 + A\] 



(12) 



[E 2 - m 2 }R{r) = 



.32 _ l d + (7 2 + ^7) _ 2i(mA 2 + A l E) 



R{r). 



(13) 



With the substitution R(r) = r 1 f/(r), it reads 



[E 2 - m 2 ]U{r) 



2 (72+^7) 

r ' o 



2i(mA 2 + AiE) 



U{r) 



(14) 



Evidently this equation is nearly identical to that of the non-Hermitian and VT —symmetric 
radial Schrodinger - Coulombic one. Of course, with the irrational angular momentum quan- 
tum number £ = — l/2 + 7 + c^/2 > 0. Its solution can therefore be inferred from the known 
non-relativistic VT— symmetric Coulomb problem ( c.f.,e.g. Mustafa and Znojil [11] and 



Znojil and Levai [15] for more details on this problem). That is 



E 2 -m< 



1/2 



n = [mA 2 + A X E] ; h = n r + 1 + 1 > 0. 



(15) 



This in turn implies 



E 

m 



AiA 2 
n 2 -A\ 



± 



A X A 2 



+ 



(n 



?,2 



Al) 



n 2 - A\ 



1/2 



(16) 



with n = n — j — 1/2 + 7, where n r = n — £+ 1 is the radial quantum number, n the principle 
quantum number, and £ = j + u/2 is the angular momentum quantum number. 

In connection with the result in equation (16), several especial cases should be inter- 
esting for they reveal the consequences of the above complexified non - Hermitian Dirac 
Hamiltonian: 

• Case 1: For A 2 = 0, the complexified Coulomb energy V(r) = —iA 1 /r = —iZa/r 
(a ~ 1/137) represents, say, the interaction energy of a point nucleus with 
an imaginary charge iZ e and a particle of charge — e. In this case 7 = 



j + 1/2) 2 + (Za) 2 , and 



E 

— = + 

m 



1 - 



{Zaf 
(n-j -I/2 + 7)' 



-1/2 



(17) 



where the negative sign is excluded because negative energies would not fulfill 
equation (15). For a vanishing potential ( Z — ) the energy eigenvalue is m. 
Obviously, unlike the ordinary ( Hermitian) Sommerfeld fine structure formula, 
equation (17) suggests that a continuous increase of the coupling strength Za from 
zero pushes up the electron states into the positive energy continuum, avoiding 
herby the energy gap. Nevertheless, for states with n — j + 1/2 one obtains 
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m V w 

The ratio E/m in (18) is plotted in figure 1 for n — 1, 2, 3, .., 10, 20, ..50. It is 
evident that as n — > oo the ratio E/m — > 1. 



Case 2 : For Ai = 0, 7 = V-^ 2 — ^b 2 and equation (16) reads 



* = ± 

m 



1 + ^ 



(n - j - 1/2 + 7) 5 



1/2 

(19) 



In this case both signs are admissible and thus two branches of solutions exist, 
but not in the energy gap. The solutions of positive and negative energies ex- 
hibit identical behaviour, which reflects the fact that scalar interactions do not 
distinguish between positive and negative charges. Moreover, states with nega- 
tive energies are pulled down to dive into the negative energy continuum, while 
states with positive energies are pushed up to dive into the positive energy con- 
tinuum. Yet, the flown away states phenomenon reemerges and for A 2 = \K\ 
states with n = j + 1/2 fly away and disappear from the spectrum. Of course one 
should worry about the critical values of the coupling (i.e., A 2 , cr it = \K\) where 
imaginary energies would be manifested. 

Case 3 : For A 1 = A 2 = A, 7 = \K\ and 

E A 2 h 2 

= ~„ 777 ± 777, 20 

m n 2 - A 2 n 2 - A 2 

Obviously, the negative sign must be discarded for it implies E = —m and thus 
contradicts equation (15). Hence, equation (20) reduces to 

— 1 + 4^- (2D 
m n 2 - A 2 

Part of this spectrum (i.e. for the principle quantum number n = 1,2,3, ..,6 ) 
is plotted in figure 2. As the coupling strength A increases from zero to n, the 
electron states are pushed up from E = m into the positive energy continuum 



avoiding the energy gap between — m to m. However, all states with n = A fly 
away and disappear from the spectrum. Nevertheless, as A increases from n and 
at A — > oo all energy states cluster just below E = —m. 

Case 4 ■' If w e replace (7 2 + £7) with 1(1 + 1), where £ = —1/2 + 



(£+ 1/2) 2 + A\ — A\ , equation (14) reduces to Klein-Gordon [23] with com- 
plex Coulomb-like Lorentz scalar and Lorentz vector potentials, S(r) = —iA 2 /r 
and V(r) = —iAi/r, respectively. That is 



[E 2 - m 2 ]U(r) 



;2 £(£+1) 2i (mA 2 + AiE) 



-di + 



U{r). 



(22) 



Which when compared with the non-Hermitian VT — symmetric Schrodinger- 
Coulomb equation implies that 



E 2 — m 2 



V2 ~ 



N = [mA 2 + AiE] ; N = n r + £ + 1 > 0, 



(23) 



and 



E A X A 2 



m 



N 2 -A 2 



± 



AiA 2 V (N 2 + A 2 ) 



+ 



s N 2 -A 2 J N 2 -A 2 



1/2 



This in turn, following similar analysis as above, yields 

; N — n — £ — 1/2 + + 1/2) 2 + A 2 , 



E 

— = + 

m 



j2 -1-1/2 
l-4i 



AT2 



for A 2 = and Ai 7^ 0, 



— = ± 

m 



1 + 

N 2 



1/2 



for Ai = and A 2 7^ 0, and 



E 

— = + 



1 + 



2 A 2 



n 2 - A 2 



■ N = n, 



(24) 



(25) 



; N = n-£-l/2 + ^/(£+l/2) 2 -A 2 , (26) 



(27) 
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for Ai = A 2 = A. Clearly, spin-0 states follow similar scenarios as those for spin- 
1/2 states ( i.e., e.g., flown away states, pushed up into the positive continuum 
and/or pulled down into the negative continuum ..etc.). 

To summarize, we have used a similarity transformation to extract exact energies 
for Dirac - particle in the generalized complex Coulomb potential. Whithin such 
non-Hermitian settings we have also obtained exact energies for Klein-Gordon 
particle. 
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Figures Captions 



Figure 1. The ratio E/m of (18) at different vales of Za for the states (from top to 
bottom) with the principle quantum number n — 1, 2, 3, 10, 20, 30, 40, and 50. 

Figure 2. Part of the spectrum E/m of (21) at different values of the coupling A and 
for states with n = 1, 2, 3, .., 6. 
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